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Abstract We discuss the application of the density functional theory in the local density 
approximation (LDA) near a ferromagnetic quantum critical point. The LDA 
fails to describe the critical fluctuations in this regime. This provides a finger- 
print of a materials near ferromagnetic quantum critical points: overestimation 
of the tendency to magnetism in the local density approximation. This is in 
contrast to the typical, but not universal, tendency of the LDA to underestimate 
the tendency to magnetism in strongly Hubbard correlated materials. We pro- 
pose a method for correcting the local density calculations by including critical 
spin fluctuations. This is based on (1) Landau expansion for the free energy, 
evaluated within the LDA, (2) lowest order expansion of the RPA susceptibility 
in LDA and (3) extraction of the amplitude of the relevant (critical) fluctua- 
tions by applying the fluctuation-dissipation theorem to the difference between 
a quantum-critical system and a reference system removed from the quantum 
critical point. We illustrate some of the aspects of this by the cases of NisAl 
and NisGa, which are very similar metals on opposite sides of a ferromagnetic 
quantum critical point. LDA calculations predict that NiaGa is the more mag- 
netic system, but we find that due to differences in the band structure, fluctuation 
effects are larger in NiaGa, explaining the fact that experimentally it is the less 
magnetic of the two materials. 

Keywords: quantum criticality, magnetism, density functional theory, first-principles calcu- 
lation. 
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Introduction 



Recent low temperature experiments on clean materials near ferromagnetic 
quantum critical points (FQCP) have revealed a remarkable range of unusual 
properties, including non-Fermi liquid scalings over a large phase space, un- 
usual transport, and novel quantum ground states, particularly coexisting fer- 
romagnetism and superconductivity in some materials. Although criticality 
usually implies a certain universality, present experiments show considerable 
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material dependent aspects that are not well understood, [1] e.g. the differences 
between UGe2 and URhGe [2, 3] and ZrZn2, [4] which both show coexisting 
ferromagnetism and superconductivity but very different phase diagrams, in 
contrast to MnSi, where very clean samples show no hint of superconductivity 
around the QCP, possibly because of the lack of the inversion symmetry. [5] 

Moreover, by far not every magnetic material can be driven to a QCP by 
pressure or by other means of supressing ferromagnetism. Typically, the tran- 
sition becomes first order as the Curie temperature, Tc is depressed. If this 
happens too far away from the fluctuation dominated regime, nothing interest- 
ing is seen. Also, more pedestrian effects are often important. For example, 
impurities or other defects can lead to scattering that smears out the quantum 
critical region. 

2, The LDA Description Near a FQCP 

One of the fingerprints of a FQCP, maybe the most universal one, is a 
substantial overestimation of the tendency to magnetism in conventional den- 
sity functional theory (DPT) calculations, such as within the local density ap- 
proximation (LDA). Generally, approaches based on density functional theory 
(DFT) are successful in accounting for material dependence in cases where 
sufficiently accurate approximations exist. Density functional theory is in prin- 
ciple an exact ground state theory. It should, therefore, correctly describe the 
spin density of magnetic systems. This is usually the case in actual state of the 
art density functional calculations. However, common approximations to the 
exact density functional theory, such as the LDA, may miss important physics 
and indeed fail to describe some materials. A well know example is in strongly 
Hubbard correlated systems, where the LDA treats the correlations in an or- 
bitally averaged mean field way and often underestimates the tendency towards 
magnetism. 

Overestimates of magnetic tendencies, especially in the LDA, are consid- 
erably less common, the exceptions being materials near magnetic quantum 
critical points (QCP); here the error comes from neglect of low energy quan- 
tum spin fluctuations. In particular, the LDA is parameterized based on the 
uniform electron gas at densities typical for atoms and solids. However, the 
uniform electron gas at these densities is stiff against magnetic degrees of free- 
dom and far from magnetic QCP's. Thus, although the LDA is exact for the 
uniform electron gas, and therefore does include all fluctuation effects there, 
its description of magnetic ground states in solids and molecules is mean field 
like. This leads to problems such as the incorrect description of singlet states in 
molecules with magnetic ions as well as errors in solids when spin fluctuation 
effects beyond the mean field are important. In solids near a QCP, the result is 
an overestimate of the magnetic moments and tendency toward magnetism {i.e. 
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Table 1. Some materials near a FQCP that we have investigated by LDA calculations. Type 
1 materials are ferromagnetic both in the calculations and in experiment; magnetic moments in 
//B per formula unit are given. Type 2 are ferromagnetic only in the calculations (calculated 
moments given, and t5?pe 3 are paramagnetic (susceptibility in 10~*emu/mol is given). The 
references are to the LDA calculations. 



Material 


ZrZn2 


NisAl 


Scsin 


FeAl 


NisGa 


Type 


1 


1 


1 


2 


2 


Calc/Exp 


0.72/0.17 


0.71/0.23 


1.05/0.20 


0.80 


0.79 


Ref. 


[16] 


[17] 


[18] 


[19] 


[17] 




Material 


Sr3Ru207 


SrRhOs 


Nao.5Co02 


Pd 




Type 


2 


2 


2 


3 




Calc/Exp 


0.80 


0.9 


0.5 


11.6/6.8 




Ref. 


[8] 


[20] 


[10] 


[21] 





misplacement of the position of the critical point) due to neglect of the quantum 
critical fluctuations. [6, 7] Examples include three types of materials: param- 
agnets that are ferromagnetic in the LDA, ferromagnets where the equilibrium 
magnetic moment is substantially overestimated in the LDA, and paramagnets 
where the paramagnetic susceptibility is substantially overestimated. 

We list examples of materials in all three categories in Table 1. At least 
two of these are cases where a large deviation between the LDA and experi- 
mental magnetic properties were noted, followed by transport measurements 
that suggest a nearby ferromagnetic quantum critical point. In particular, in 
Sr3Ru207, LDA calculations with the experimental crystal structure found a 
sizeable moment, [8] while experimentally the material was known to be a 
paramagnetic metal. Grigera and co-workers then showed that Sr3Ru207 has a 
metamagnetic quantum critical point at moderate field. [9] Pd metal provides 
another example: the calculated LDA magnetic susceptibility is nearly twice 
larger that the experimental one. Correspondingly, Nicklas et al[\Y[ found a 
FQCP in the Pdi_a;Nia; system at x = 0.026, where the transport properties 
become non-Fermi liquid. 

We emphasize that substantial overestimates of the tendency of metals to- 
wards ferromagnetism within the LDA is a rare occurance, and propose that 
it be used as an indicator of critical fluctuations in a material. However, for 
this to be an effective screen, competing states, like antiferromagnetism need 
to be ruled out in each material. An interesting case study is LiV204, which 
is a paramagnetic metal and occurs in the cubic spinel structure. Remarkably, 
it was discovered by Kondo and co-workers that this material behaves at low 
temperature like a heavy fermion metal. [12] LDA calculations showed that 
the material is unstable against ferromagnetism with a sizeable moment. [14, 
15, 13] But calculations also show that the interactions are antiferromagnetic. 
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and as a result it is more unstable against antiferromagnetism, which however 
is frustrated on the spinel lattice. While LiV204 may be near an antiferromag- 
netic QCP, it is not a material near an FQCP. 

3. "Beyond-LDA" Critical Fluctuations 

A popular way to add quantum or termal fluctuation to a mean-field type 
theory is via fluctuation corrections to Ginzburg-Landau expansion of the free 
energy. For a detailed discussion we refer the reader to the book of Moriya 
[22] and the review article of Shimizu [23]. In short, one writes the free energy 
(or the magnetic field) as a function of the ferromagnetic magnetization, M, 

EstaticiM) = ao + ^^asnM^", (1) 

n>l 

HstaticiM) = 5^a2„M2"-i (2) 

n>l 

(obviously, 02 gives the inverse spin susceptibility without fluctuations), and 
then assume Gaussian zero-point fluctuations of an r.m.s. magnitude ^ for 
each of the d components of the magnetic moment (for a 3D isotropic mate- 
rial like Pd, d = 3). After averaging over the spin fluctuations, one obtains a 
fluctuation-corrected functional. The general expression can be written in the 
following compact form: 



H{M) = ^d2nM2"-l 



n>l 



«2n 

For instance, 



i>0 



- ^ >-2 a4 a6 

0.2 = a2 + -04? + Y«6C + y"8? ■•• 

04 = 04 yOeC^ + 2l08?^... 



(4) 



The unrenormalized coefficients can be taken from fixed spin momen LDA 
calculations, in which case ? becomes the amphtude of those fluctuations only, 
which are not taken into account in LDA (as mentioned, LDA includes some 
quantum fluctuation, specifically short-range fluctuations present in the in- 
teracting uniform electron gas). In principle, one can estimate ? from the 
fluctuation-dissipation theorem, which states that (see, e.g., Refs. [24, 25]) 
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where x(q) ^) is the magnetic susceptibility and Q. is the Brillouin zone vol- 
ume. It is customary to approximate x(q, uj) by its small q, small oj expansion 
[24, 25]: 

Xo(q,w) = N{EF)-aq^ +ihu/q (6) 

With the expansion (6) the integrations can be performed analytically, and 
the final result reads: 

e = ^-^^^f^ [Q' + Q-') + HI + Q% (8) 



where Q = qc^/a/bvp, and qc is the cutoff parameter for momentum integra- 
tion in Eqn. 5 (the frequency integration at a given q is usually assumed to be 
cut off at 6<; = vpq)- 

To proceed along these lines one needs to find a way to calculate the crucial 
parameters of the expansion (6). It was suggested by Moriya [22] that these can 
be expressed as certain integrals over the Fermi surface, by expanding the RPA 
expression for xo- Below, we offer a derivation equivalent to that of Moriya, 
but rendering the results in more computable form. We start with the RPA 
expressions for the real and imaginary parts of xo '■ 

Rexo(q,0) = 5] [/(i^k) - /(£^k+q)] (^^k+q - -Bk)-' (9) 

k 

Imxo(q,^) = ^[/(£;k)-/(^k+q)]<^(^k+q-^k-a;), (10) 

k 

where f{E) is the Fermi function, —^^j§^ = 5{E — Ep). Expanding Eqn. 9 
in A = i?k+q - -Sk = Vk-q+5 Y^ap l^k^QaQp + we get to second order in 

q 

+ 6 [ dEj. ) 

The odd powers of Vk cancel out and we get (a, /3 = x, y, z) 
Rexo(q) = 

N{Ef) + 1^-^ + E-6 ^ 

a,l3 a, 13 ^ 

q^ d{N{EF)^i,,) q^ d?{N{EF)vl) 

= ^^^^^ + T — dE; — + ¥ — dEi — ' 
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where v"^ = Vy = v^, = fJ^yy = l^zz- The last equality assumes cubic 
symmetry; generalization to a lower symmetry is trivial. Using the following 
relation, 



one can prove that 



(f{N{EF)vl)_ d{N{EF)lixx) 



dEl 



Therefore 



Rexo(q) =iV(^F)- 



dEp 

(P {N{Ef)vI) 
12 dEl 



Similarly, for Eqn. 10 one has 

Imxo(q,'^) = E 

k 



df{e) 
de 



a;^(vk-q — oj) 



(11) 



(12) 



(13) 



After averaging over the directions of q, this becomes, for small uo, 



ImXo(g,t^) 



9 E -r-r^(^kQ' 



u;) = ^{NiEF)v- 



(14) 



Although in real materials the Fermi velocity is obviously different along dif- 
ferent directions, it is still a reasonable approximation to introduce an aver- 
age VF- Then the above formulae reduce all parameters needed for estimat- 
ing the r.m.s. amplitude of the spin fluctuations to four integrals over the 
Fermi surface, specifically, the density of states, N{Ef), a = {^i^F)vJ ^ 



,(N{Ep)vl) 



b=\ {N{Ef)v 1) and - y ^^nEpY 

The physical meaning of these parameters is as follows, a defines the rate 
at which the static susceptibiUty x{Qi 0) f^Us away from the zone center, i.e. 
the extent to which the tendency to ferromagnetism is stronger than that to 
antiferromagnetism. This translates into the phase space in the Brillouin zone 
where the spin fluctuations are important, h controls the dynamic effects in 
spin susceptibility. 

Note that the cutoff parameter Qc remains the only undefined quantity in this 
formalism. One obvious choice is Qc = ^N{EF)/a, because for larger q the 
approximation (6) gives unphysical negative values for the static susceptibiUty. 
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On the other hand, one may argue that Qc should reflect mainly the geometry 
of the Fermi surface and thus not depend on a at all. We will come back to 
this issue later in this paper and will propose an approach that avoids using qc 
whatsoever. 

4. NisAI and NiaGa 

Here we use the closely related compounds NisAl and NisGa to illustrate 
some of the above ideas. Further details may be found in Ref. [17]. These 
have the ideal cubic CuaAu cP4 structure, with very similar lattice constants, 
a = 3.568 A and a = 3.576 A, respectively, and have been extensively studied 
by various experimental techniques. NisAl is a weak itinerant ferromagnet, 
Tc = 41.5 K and magnetization, M=0.23 /is/ceU (0.077 /xs/Ni atom) [26] 
with a QCP under pressure at Pc=S.l GPa, [27] while NisGa is a strongly 
renormaUzed paramagnet. [28] Further, it was recently reported that NisAl 
shows non-Fermi liquid transport over a large range of P and T range down to 
very low T. [29] 

Previous LDA calculations showed that the magnetic tendency of both ma- 
terials is overestimated within the LDA, and that NisGa is incorrectly predicted 
to be a ferromagnet. [30-35] Moreover, in the LDA the tendency to magnetism 
is stronger in NisGa than NisAl, opposite to the experimental trend. This poses 
an additional challenge to any theory striving to describe the material depen- 
dent aspects of quantum criticality. The two materials are expected to be very 
similar electronically (the small difference between the two is due to relativis- 
tic effects associated with Ga in NisGa). Thus these two very similar metals 
offer a very useful and sensitive benchmark for theoretical approaches. We use 
this to test an approach based on the fluctuation dissipation theorem applied 
to the LDA band structures with an ansatz for the cut-off Qc. We find that this 
approach corrects the ordering of the magnetic tendencies of the materials, and 
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gives the right ground states at ambient pressure as well as a reasonable value 
of Pc for NigAI. 

The LDA calculations were done using the general potential linearized aug- 
mented planewave (LAPW) method with local orbital extensions [36, 37, 39] 
as decribed in Ref. [17], with the exchange-correlation functional of Hedin and 
Lundqvist with the von Barth-Hedin spin scaling [40, 41]. The LDA electronic 
structure is given in Fig. 1 and Table 2, while results of fixed spin moment cal- 
culations of the magnetic properties at the experimental lattice parameters and 
under hydrostatic compression are given in Figs. 2 and 3. The two compounds 
are very similar in both electronic and magnetic properties, the main apparent 
difference being the higher equilibrium moment of NiaGa (0.79 /is/f.u. v*. 
0.71 )UB/f.u.), in agreement with other full potential calculations. [34, 35] 

The propensity towards magnetism may be described in terms of the Stoner 
criterion, IN{Ep), where / is the so-called Stoner parameter, which derives 
from Hund's rule coupling on the atoms. For finite magnetizations, the so- 
called extended Stoner model [42], states that, to the second order in the spin 
density, the magnetic stabilization energy is given by 

r-M 

AE = M^[ mdm/2N{m) - 1/4], (15) 
Jo 

where N{M) is the density of states averaged over the exchange splitting cor- 
responding to the magnetization M. Fitting the fixed spin moment results to 
this expression, we find Iai = 0.385 eV and loa = 0.363 eV. These gives 
IN{Ef) =1.21 and IN{Ef) = 1.25 for NisAl and NiaGa, respectively. Both 
numbers are larger than unity, corresponding to a ferromagnetic instability, and 
the value for NiaGa is larger than that for NiaAl. Importantly, the difference 
comes from the density of states, since Iai > ^Ga- In both compounds, mag- 
netism is suppressed by compression, with an LDA critical point at a value 
5a/ a ^ -0.05 - -0.06. In NiaAl, the critical point at 5a /a =-0.058 corresponds 
to the pressure of Pc =50 GPa, [43] which is much higher than the experimen- 
tal value. It is interesting that, as in ZrZn2 [16], the exchange splitting is very 
strongly k-dependent; for instance, in NiaAl at some points it is as small as 40 
meV/fi B near the Fermi level, while at the others (of pure Ni d character) it is 
close to 220 meV///B. 

Notwithstanding the general similarity of the two compounds, there is one 
important difference near the Fermi level, specifically, the light band crossing 
the Fermi level in the middle of the F-M or F-X directions is steeper in NiaAl 
(Fig. 1). This, in turn, leads to smaller density of states. This comes from 
a different position of the top of this band at the V point, 0.56 eV in NiaGa 
and 0.85 eV in NiaAl. The corresponding electronic state is a mixture of Ni p 
and Al (Ga) p states, and is the only state near the Fermi level with substantial 
Al (Ga) content. Due to relativistic effects, the Ga p level is lower than the 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 2. Energy vs. fixed spin moment for NisAl and NisGa at the experimental lattice 
parameters. The energy zero is set to the non-spin-polarized value. 

Table 2. Magnetic energy (see text), magnetic moment in hb/c&U and N{Ef) in eV~^ for 
NisAl and NisGa on a per spin per formula unit basis. 

\AE\ (meV) M (calc.) M(expt.) N{Ef) 

NisAl 103 OTI 023 3.2 

NisGa 14.3 0.79 0.00 3.4 

Al p level and this leads to the difference in the position of the corresponding 
hybridized state. Note that this is a purely scalar relativistic effect. Including 
spin orbit does not produce any further discernible difference. 

Returning to magnetism, the fixed spin moment calculations provide the 
energy as a function of the magnetization M (Fig. 2). One can write a 
Landau expansion for E[M) as in Eqn. 1, which may then be treated as a 
mean field expression adding the effects of spin fluctuations. [23] 

Treating this as a mean field expression and adding the effects of spin fluc- 
tuations [23] leads to renormalization of the expansion coefficients. The renor- 
malized coefficients are written as power series in the averaged square of 
the magnetic moment fluctuations beyond the LDA, as in Eqn. 3. ^ may 
then be estimated by requiring that the corrected Landau functional reproduces 
the experimental magnetic moment (for NiaAl) or experimental magnetic sus- 
ceptibility (for NisGa). The "experimental" ^'s obtained in this manner are are 
0.47 and 0.55, respectively, which imphes that spin fluctuation effects must be 
stronger in NisGa than in NisAl. 

A link can now be made between this fact and the electronic structures, 
using the formalism outlined in the previous section. As discussed, the cutoff 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Figure 3. FSM calculations under hydrostatic pressures. Magnetic energy, defined as the 
energy relative to the non- spin-polarized result at the same volume, as a function of the moment 
and linear compression. Left and right panels correspond to NisAl and NisGa, respectively. 



parameter Qc is the least well defined quantity in this formalism. Furthermore, 

the fermiology of these compounds is very complicated: in the paramagnetic 
state, there are four Fermi surfaces, two small and two large (one open and one 
closed). In this situation, it is hardly possible to justify any simple prescription 
for Qc. Therefore, we chose a different route: we assume that Qc is the same 
for both materials, and choose a number which yields a good description of 
both the equilibrium moment in NiaAl and the paramagnetic susceptibility in 
NisGa, Qc = 0.382 Oq ^. Note that this is larger that the diameters of the small 
Fermi surfaces but smaller than the radius of the Brillouin zone, » 0.5 Aq ^. 

To calculate the above quantities, especially a, we need accurate values of 
the velocities on a fine mesh. Numerical differentiation of energies within the 
tetrahedron method proved to be too noisy. Therefore we use the velocities 
obtained analytically as matrix elements of the momentum operator, computed 
within the optic program of the WIEN package. A bootstrap method, [44] as 
described in Ref. [21], was used to obtain stable values for a, b. We found 
for NisAl, using as the energy unit Ry, the length unit Bohr, and the velocity 
unit Ry Bohr, a = 230, b = 210, vp = 0.20, and ^ = 0.445 hb- For NiaGa 
a = 140, b = 270, vp = 0.19, and ^ = 0.556 hb- Using the resulting values 
of ^ each compound we obtain a magnetic moment of M = 0.3 hb/ccW for 
NiaAl and a paramagnetic state with the renormalized susceptibility x(0, 0) = 
l/a2 = 6.8 X 10~^ emu/g for NisGa, thus correcting the incorrect ordering of 
the magnetic tendencies of these two compounds and reproducing extremely 
well the experimental numbers of M = 0.23 hb, x(0) 0) = 6.7 x 10^^ emu/g, 
respectively. This qualitative behavior is due to the different coefficient a, i.e., 
different q dependencies of xo{q, 0) at small q, which relates to the phase space 
available for soft fluctuations. 



Density Functional Calculations near Ferromagnetic Quantum Critical Points 1 1 



Now we turn to the pressure dependence. The above results imply that 
beyond-LDA fluctuations are already larger than the moments themselves at 
P = 0. In this regime, we may assume that the size of the beyond-LDA fluc- 
tuations is only weakly pressure dependent. Then we can apply Eqn. 3 to the 
data shown in Fig. 3 using ^ = 0.47 as needed to match the P = value of M. 
This yields a value Pc=10 GPa in quite good agreement with the experimental 
value, Pc=8.1 GPa. [27] 

5. Towards a Fully First Principles Theory 

The results for NisAl and NisGa, discussed above, and in Ref. [17], show 
that an approach based on correction of the LDA using the fluctuation dissipa- 
tion theorem has promise. However, the results hinge on an unknown cut-off, 
which serves the purpose of including fluctuations that are associated with the 
FQCP and are not included in the LDA, from those that are included in the 
LDA. While it is apparently possible to obtain useful results using reasonable 
ansatz for this cut-off, it would be much better to have a truly first principles 
theory, with no parameters. In order to construct such a theory, one should find 
a way of solving the double counting problem, i.e including in the correction 
only those fluctuations that are not already taken into acount at the LDA level. 
This amounts to subtracting from Eqn. 5 the fluctuations already included 
in the LDA. Since the LDA is known to work well for materials far from an 
FQCP, this means that the correction should be zero or close to it for the most 
materials. 

We suggest that a consistent way to accomplish this is by introducing a 
"reference" susceptibiUty x^^'^(Q)'*^) subtracting it from x(q,a;) : 

e = ^ jd'q |^llm[x(q,^) -x'-^^lq,^)], (16) 

We shall use the same expansion6 for both x(q, uj) and x^'^H'i' '^)' ^ derive 
equivalent expansions 

X-'(q, iv) = XoHO, 0) - / + V - iBu/q, (17) 

where Xo^(0, 0) = l/N{Ep) (density of states per spin) is the bare (non- 
interacting) static uniform susceptibiUty, and the Stoner parameter / is only 
weakly dependent on q and u. Note that A = a/N'^, B = h/N"^ , where a 
and b are the coefficients introduced in Eq.6. We also introduce a notation, 
A = N{Ep)~^ — I. As long the same functional form (17) is used for x(q, 
and x^^'^ i^.:^), the condition for the convergence of the integral (16) is that 
the coefficients A and B, controlling the short-range and high frequency fluc- 
tuations are the same. Of course, the parameter A, defining the proximity to 
the QCP, is different in the reference system, which like the uniform electron 
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gas upon which the LDA is based, should be far from any QCP (let us call A 
for the reference system Aq). 

To calculate the integral ((16), we write it in the following form: 

= f /^''^/ ^^Mx(A,q,u;)-x(Ao,q,u;)]. (18) 

For instance, x(0,q, w) represents the susceptibility right at the FQCP. This 
diverges for q = 0, w = 0. The derivation then proceeds as follows: 



Where we introduce the Landau cutoff frequency, ooc = vq (here v is an aver- 
age Fermi velocity) and the notation /3 = Bv. We will also need the following 
function: 

F(A,/3,x) = j a^'dx\n[{^ + x^f + 0^] 
(A + x2)2+/32 



4 



^i^±^{ln[(A + x' f + - 2} + /?Atan-^ 



Now 



e = ^lunjFiA,P,Q)-FiA,,p,Q) 

-F(A, 0, Q) + F(Ao, 0, Q) - F{A, /3, 0) 
+F(Ao, /3, 0) + F(A, 0, 0) - F(Ao, 0, 0)]. 

This is particularly easy to evaluate at A = 0. The result is 

S^(Ao) = 

2 ,,7r lAo, /3^-A2 A^ + f]^ A^ 

-[Ao/3(- - tan-i -^) + ^—-—^ In ° / + In 



^\So) = -^^[ASote.n-\So') + In(l + 5^) - ^^Inll + ^o'^)], 

where = AoN^/bvp. Obviously, for arbitrary A the answer is simply 

e = E\So)~E\S). (20) 



Given that usually the reason for a quantum criticality is a large density of 
states, it makes sense to take the Stoner parameter for the reference system 
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the same as for the system in question. The point is that the density of states 
is a highly non-local parameter (note that it involves a delta function integral 
in energy), which can hardly be discerned from local information about the 
charge density, while the Stoner parameter is a very local quantity associated 
with the exchange-correlation potential. The difference between A and Aq 
then comes from the difference between N = N{Ep) and the density of states, 
A^O) of the reference system. 

One may think about several different ways for choosing A^o- One may be 
to take average N{E) over the width of the valence band, A^o = n/t, where n 
is the total number of states in the band and t is its width. One can also think 
about the density of states of the uniform electron gas with the same Stoner 
parameter. There may be other, more sophisticated prescriptions. Probably, 
the most practical approach will be found after several trial and error tests with 
real materials. 

6. Summary and Open Questions 

The failure of the usual approximations to density functional theory, for 
example, the LDA, to describe the magnetic properties of materials near fer- 
romagnetic quantum critical points is associated with renormalization due to 
critical fluctuations. It is pointed out that since such fluctuations are invari- 
ably antagonistic to ferromagnetic ordering, deviations between experiment 
and LDA calculations in which the LDA is overly ferromagnetic can be a use- 
ful screen for materials near FQCPs. These errors in the LDA can be corrected 
using a phenomenalogical Landau function approach with the fluctuation am- 
plitude as a parameter. However, there is hope that this parameter can be ob- 
tained from the electronic structure via the fluctuation dissipation theorem and 
a suitable reference system. The key remaining challenges in our view are to 
define the reference system to be used, and to use calculations to determine the 
usefulness of this approach for real materials near a critical point. 
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